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ABSTRACT 

Three dimensional (3D) adaptive-mesh refinement (AMR) hydrodynamical simulations of the wind- 
wind collision between the enigmatic super-massive star 77 Car and its mysterious companion star 
are presented which include radiative driving of the stellar winds, gravity, optically- thin radiative 
cooling, and orbital motion. Simulations with static stars with a periastron passage separation reveal 
that the preshock companion star's wind speed is sufficiently reduced that radiative cooling in the 
postshock gas becomes important, permitting the runaway growth of non-linear thin shell (NTSI) 
instabilities which massively distort the WCR. However, large-scale simulations which include the 
orbital motion of the stars, show that orbital motion reduces the impact of radiative inhibition, and 
thus increases the acquired preshock velocities. As such, the postshock gas temperature and cooling 
time see a commensurate increase, and sufficient gas pressure is preserved to stabilize the WCR 
against catastrophic instability growth. We then compute synthetic X-ray spectra and lightcurves 
and find that, compared to previous models, the X-ray spectra agree much better with XMM-Newton 
observations just prior to periastron. The narrow width of the 2009 X-ray minimum can also be 
reproduced. However, the models fail to reproduce the extended X-ray mimimum from previous 
cycles. We conclude that the key to explaining the extended X-ray minimum is the rate of cooling of 
the companion star's postshock wind. If cooling is rapid then powerful NTSIs will heavily disrupt the 
WCR. Radiative inhibition of the companion star's preshock wind, albeit with a stronger radiation- 
wind coupling than explored in this work, could be an effective trigger. 

Subject headings: hydrodynamics - stars:early-type - X-rays:stars - stars:binaries - stars:winds, outffows 
- stars individual (7^ Carinae) 



1. INTRODUCTION 

Of the known massive stars in our galaxy, 77 Car is pos- 
sibly the largest and finest example of a pre-hypernova 
candidate, presenting a rare but exceptional opportu- 
nity to test our current understanding of stellar evolu- 
tion in the upper Hertzsprung- Russell diagram. Yet, 
this gem is not without its fiaws as, since the "Great 
Eruption" which formed the Homunculus nebula and the 
mor e recent eruption wh ich formed the Little Homuncu- 
lus (jlshibashi et al.|[2QQ3l ). r] Car has been enshrouded by 
a dusty cocoon which complicates observations in the UV 
and optical. Fortunately, r] Car is extremely bright at X- 
ray wavelengths which suffer less from extinction, thus 
providing an invaluable probe of the inner nebula. 
The except ional monitoring of 7 7 Car at X-ray 
^ 3] IlMll, [2003: ICorcoranI 



19991: iPittard fc CorcoranI 120021: lAkashi et aP I20Q6I: 
Henleve^f] 120081: lOkazaki et all i2008i : iParkin et al.l 



wave le ngths ([Corcoran et . 

20051: iCorcoran et all [2Q1Q|; Eamaguchi et aP 120071: 



Henlev et all 120081 : iLevder et alJIlooM 120101 : iPian et al.l 
20091 ) has characterised the periodic variability as in- 
dicative of a highly eccentric (e 0.9), long-period 
(~ 5.54 yr) binary system (Table [ij). The bi- 
nary hypothesis a ppears to be well supported by 
theoretical models (jPittard et al.l 119981 : llshibashi et al.l 

|parkin@mso.anu.edu.aul 



2009: iK^i Sokei]l2009bD. as we ll as by observations 
at infrared ("Whitelock et al.' "1994", "2004"; "Nielsen et_al| 
2009), radio (Duncan & White 2003; Abraham e tHI 
2 005 1)1). optical (iDaminelii 1996; Damineli et a l.l I200QL 
l2008J bl: Ivan Gende ren et al. 2003; Nielsen et alTl2007[ ). 
far ultra -violet (Ep ing et al. 2005), and 7— ray wave- 
lengths (iTavani et al., .2009.: .Walter et al.l I2010L - al- 
though seelOhmeiaD Eoifl and^^ In 
this scenario the X-ray emission originates from the shock 
heated plasma generated by the fast wind of a companion 
star which ploughs into the slow dense wind of 7^CarEl 
(e.g. IPittard et aIlll998D . 

The observed X-ray lightcurve displays the general 
characteristics expected from a high orbital eccentric- 
ity colliding winds binary, but has some puzzling fea- 
tures. For instance, the most recent X-ray minimum was 
considerably short er than the extended minima observed 
in previous cycles ([Corcoran et al.|[2001l : lCorcoranll2005l : 

^ For the remainder of this paper we refer to the larger star 
(commonly referred to as a luminous blue variable (LBV) star - 
although see Mehner et al. 2010a) and the smaller comp anion star 
(es timated to be of early O type or a WR star - P ittard &: CorcoranI 
[2002: Verner et al. 2005; Teodoro et al. 2008; Mehner et al.ll2010al ) 
as the primary and companion star, respectively. 
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iCorcoran et alJ[2QlQf ) . The exact cause of the X-ray min- 
imum is a matter of debate (for an o verview of some of 
the possibilities see lParkin et aljr2QQ9f ). yet it is now clear 
that models must also be able to account for the observed 
cycle-to-cycle variation. 

The nature of the wind-wind collision region (WCR) 
has recently been investigated with 3D models. Us- 
ing a Smoothed-Particle Hydrodynamics (SPH) model 
of the wind-wind collision, Okazaki et al. (2008) found 
that they could match the extended X-ray minimum. 
In contrast, adopting t he 3D dynamical m odel of 
iParkin fc PittardI (j2008[ ). iParkin et all (j2009[ ) showed 
that when the spatial extent of the X-ray emission region 
and energy dependence of the emission and absorption 
were taken into consideration, the width of the observed 
X-ray minimum could not be reproduced by an eclipse 
of the X-ray emitting plasma alone. Furthermore, mod- 
els with the preshock stellar winds at terminal velocity 
over-predicted the observed X-ray emission in the 7-10 
keV band by an order of magnitude, indicating that a 
reduction in the preshock speed of the companion star's 
wind is required. 

Such a reduction is plausible, as the highly eccentric 
orbit causes the WCR to enter into the wind accelera- 
tion region of the companion star around periastron pas- 
sage. Considering the immense luminosity of the primary 
star, the acceleration of the companion's wind may also 
be significantly inhibited. A reduction in the preshock 
velocity will cause radiative cooling to become increas- 
ingly important in the postshock compa nion's wind, and 
will affect the sta bility of the WCR (jDavidsonI 120021 : 
IParkin et al.ll200"9l ). In fact, a ram pressure balance be- 
tween the stellar winds may be lost and the shoc k may 
collapse onto the companion star (jParkin et al.l [2Q09f ). 
In this case the 7-10 keV flux which originates predom- 
inantly from the apex of the WCR would be quenched, 
a feature which is necessary in models aiming to explain 
the extended X-ray minimum. However, it is unclear 
why a collapse of the WCR (if one occured) would be of 
a much shorter duration in the most recent periastron 
passage whilst being extended in previous cycles. 

In this paper we describe our results from a suite of 
three dimensional (3D) simulations in which we have im- 
plemented radiative driving, gravity, orbital motion, and 
optically thin radiative cooling, allowing the role of wind 
acceleration, interacting radiation fields, and instabilities 
in the wind-wind collision region (WCR) on the gas dy- 
namics and resulting X-ray emission to be explored. We 
first present simulations with static stars at a separa- 
tion corresponding to periastron passage. The turbulent, 
highly unstable nature of the wind-wind collision region 
(WCR) is revealed, and results show that when the ac- 
celeration regions of the stars are considered the condi- 
tions in the postshock gas permit the growth of power- 
ful n on-linear thin-shell instabilities (NTSI - Vishni^ 
[1994). Subsequent vigorous oscillations lead to the col- 
lision of dense fragments of the WCR against the com- 
panion star. Importantly, the 7-10 keV X-ray luminosity 
shows a marked reduction when compared to an equiv- 
alent simulation which neglects the wind acceleration. 
We then perform large-scale, high resolution simulations 
which include the orbital motion of the stars, and find 
that the rapid orbital motion of the stars around perias- 
tron acts to: reduce the degree of radiative inhibition of 



TABLE 1 

Adopted system parameters for 77 Car. 



Parameter 



Value Reference 



Orbital period (d) 


2024 


1 


Eccentricity (e) 


0.9 


1 


a (au) 


16.64 


2 


cisep((/> = 0.0) (1013 cm) 


2.48 




dsep{(p = 0.5) (1013 cm) 


47.3 




Distance (kpc) 


2.3 


3 


ISM + nebula column (lO^^cm-^) 


5 


4 



Note. — a is the semi-major axis of the orbit and ctsep {(f)) is the 
binary separ ation as a funct i on of o rbital phase. References are as 
follows: 1 ^ 'Damineh et al.l (l2008bl V 2 =IHillier et al.l (»2001i) , 3 = 
lDavidson°^ Humphrevsl (|1997l ), T = IHamaguchi et all (|2007I V 



the preshock companion star's wind. Consequently, the 
catastrophic disruption of the WCR seen in the static 
stars simulation is not reproduced. 

We conclude that the behaviour of the observed X-ray 
emission from rj Car around periastron is tied to the rate 
of radiative cooling in the postshock companion's wind. 
The suppression of the preshock companion's wind by ra- 
diative inhibition (albeit with a stronger radiation-wind 
coupling than explored in this work) could be an im- 
portant trigger for rapid cooling. Furthermore, the dif- 
ferences between the observed short and long duration 
X-ray minima may be due to the companion's postshock 
gas conditions lying close to the dividing line between 
instability growth which merely perturbs the WCR, and 
far more vigorous NTSIs which destroy it (possibly caus- 
ing the WCR to collapse onto the companion star in the 
process). The latter case may explain the longer observed 
minima. 

The remainder of this paper is structured as follows: 
a description of the hydrodynamic model and the X-ray 
emission calculations are given in § [2l In § [3] we review 
some theoretical background relevant to our current in- 
vestigation. The results of the hydrodynamic simulations 
are presented in § [H (static stars at periastron separa- 
tion), and § [5] (large-scale orbit). A discussion and sug- 
gestions for future directions are given in § [6l and we 
close with our conclusions in § [71 

2. THE MODEL 

2.1. Hydrodynamic modelling 

The wind-wind collision is modelled by numerically 
solving the time-dependent equations of Eulerian hydro- 
dynamics in a 3D Cartesian coordinate system. The rel- 
evant equations for mass, momentum, and energy con- 
servation are: 



^+V-/5VV + VP = /9f, 

ot 

at KrriYi 



(1) 
(2) 

A(T) + pf • V. (3) 



Here E = e + ||v| , is the total gas energy, e is the 
specific internal energy, v is the gas velocity, p is the 
mass density, P is the pressure, T is the temperature, 
and mn is the mass of hydrogen, f is the force per unit 
mass and includes gravity and radiative driving terms. 
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We use the ideal gas equation of state, P = (7 — l)pe, 
where the adiabatic index 7 = 5/3. 

The radiative cooling term, A(T), is calculated from . 
the MEKAL thermal plasma code (jMewe et all 119951 : 
iKaastral [19921 ) distributed in XSPEC (vll.2.0). The 
temperature of the unshocked winds is assumed to be 
maintained at ~ 10^ K via photoionization heating by 
the stars. Throughout this work solar abundances are 
assumed (Anders & Grevesse 1989). 

The body forces acting on each hydrodynamic cell 
are the vector summation of gravitational forces from 
each star, and continuum and line driving forces from 
the stellar radiation fields. The computation of the 
line acceleration is based on a local Sobolev (1960) 
treat ment of the line transport, fo llowing the stan- 
dard ICastor. Abbott, Klein fl975)l (hereafter iCAK) 
formalism developed for single OB star winds. To ac- 
count for the correction to the line force due to the fi- 
nite size of the stellar disk (i.e. the finite disk correction 
factor - LC astor 1974; Pa uldrach et al.l I1986D we follow 
iCranmer fc Owocki (1995) and define a "wind centred" 
coordinate system {x'^y'^z') to calculate the local ve- 
locity gradient (see Fig. [T]) and then move to the "star 
centred" coordinate system {x^y^z) to determine the pro- 
jected velocity gradient along the direction vector n using 
the following transformation, 

X = {x' cos ^0 + sin ^0) cos — y' sin 0o 

y =y' cos + ix' cos 6^0 + z' sin ^0) sin (4) 

z = z' COS 6^0 — x' sin ^0 

The vector radiative force per unit mass, grad^ is then 
computed from a numerical Gaussian integration (using 8 
directional vectors) of the intensity /(n) times n- V(n-'y) 
over the solid angle covering the stellar disk. 



^rad 



/(n) 



n • V(n • v) 



(5) 



a and k are the standard ICAKI parameters, <Je is the 
specific electron opacity due to Thomson scattering, and 
Vth. is a fiducial thermal velocity calculated for hydrogen. 
Shadowing by the companion star is accounted for in our 
calculations, and in such cases only the visible part of 
the stellar disk contributes to the radiative driving force. 
The line driving is set to zero in cells with temperatures 
above 10^ K, since this plasma is mostly ionized. Further 
details concerning the line force calculations can be fou nd 
in Cranmer & Owocki (1995) and Gay ley et al.] (|1997[ ). 

The stellar winds are initiated in the instantaneously 
accelerated and radiatively driven stellar winds simula- 
tions in two slightly different ways. In the former, appro- 
priate hydrodynamic variables (i.e. p, P, v) are mapped 
into cells residing within a radial distance i^map which is 
a function of stellar separation (see Appendix [A]) . Thus, 



Mi 



47rr?|vc 



(6) 



where ri is the distance of the cell from the respective star 
and the subscript i refers to the wind being considered. 
In contrast, to initiate radiatively driven stellar winds 
we map hydrodynamic variables determined from single 
star wind profiles (calculated with the parameters noted 
in Table [2]) within a stellar radius oi c^l.l after every 



y 



Fig. 1. — Schematic diagram illustrating the transformation from 
the "wind centred" to the "star centred' coordinate system. 



time step. As noted bv iPittardI (|2QQ9[ ). it is necessary to 
resolve the region above the stellar surface (i.e 0.1 R^) 
with at least 3 cells. To distinguish between the stellar 
winds in the simulations we include an advected scalar 
variable in the hydrodynamics calculations. 

It should be noted that the primary is an LBV, and 
therefore it is unclear to what level usins: the [GAKI for- 
malism is appropriate to describe the wind driving from 
this star. Furthermore, Smit h et al.l (|2003) inferred a 
higher rate of mass-loss at the poles based on the ob- 
servation of varying absorption as a function of latitude 
in the refiected emission from the Homunculus nebula, 
which may be related to the rotation rate of the pri- 
mary star, though it m ay also be related t o the presence 
of the companion star (jGroh et al.ll2QlQ"al ). Since we do 
not consider a latitude dependence, the mass-loss rates 
which we infer are likely indicative of the equatorial wind 
(assuming the equatorial and orbital planes are aligned). 

Adopting solar abundances may also affect the i nferred 
mass-loss rates. For instance, iHillier et al.l (|2QQll ) found 
N to be significantly enhanced, while C and O could be as 
low as 1/25 and 1/50 solar, respectively. Although these 
differences will alter the opacity of the primary wind, we 
do not anticipate that the observed X-ray spectrum at 
E > 1 keV will be significantly affected. 

2.2. The hydrodynamic code 

The simulations presented in this wo rk were per- 
formed using versio n 3.1.1 of the FLASH (|Frvxell et al.l 
I2QQQI : iDubev et all I2QQ9D hydrodynamical code. This 
code operates with a block-structured AMR grid (e.g. 
iBerger OHgerl 119891) usi ng the PARAMESH pack- 
age (jMacNeice et al.l boQOl ) under the message-passing 
interface ( MPI) architecture. The pi ecewise-parabolic 
method of iColella fc Woodward (|1984f ) is used to solve 
the hydrodynamic equations. Customized units have 
been implemented into the FLASH code for radia- 
tive cooling for optica lly-thin plasma using t he ex plicit 
method described by iStrickland fc BlondinI ()1995f ). ra- 
diative driving, gravity, and orbital motion. Simulations 
are performed with the stars fixed at periastron (§[4|) and 
as they move through an orbit (§[5j), and a description of 
the grids used can be found in the relevant section (i.e. 
size, shape, resolution, and refinement criteria). 

We note that as a goal of our investigation is to ex- 
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TABLE 2 

Parameters used to calculate the line driving of the 
stellar winds. 



Parameter 



Primary Secondary Reference 



M (Mo) 
(R©) 
Tcs (K) 
L* (lO^Lo) 
k 
a 

M(M0yr-i) 
Voo ( km s"-*- ) 
Jl 



120 
100 
25,800 
4 

0.30 
0.52 

4.8 X 10- 
500 



30 
20 

30,000 
0.3 
0.50 
0.68 

1.4 X 10-5 
3000 



0.18 



Note. — M is the stellar mass, i^* is taken to be the radius of 
the gravitationally bound core of the star (which for the primary 
star is taken to be at a radius of 100 R©, i.e. the photosphere 
will exist somewhere in the stellar wind), Tcs is the temperature 
at i?* , is the stellar luminosity, M is the stellar mass-loss rate, 
Voo is the wind terminal velocity, and 77 = M2 00 2/^/^1'^ 00 1 is the 
wind momentum ratio, k and a are the CAK line driving param- 
eters, where subscripts 1 and 2 are used to define the coupling 
between the winds and the radiation fields of the primary and 
seconda ry star respectiv ely. Note that single star mCAK calcu- 
lations (|Pauldrach et al.l [l986) were performed to determine val- 
ues of k and a required to produce the desired M and Voo for 
ea ch wind. Referenc es are as fol lows: 1 = Hillier et al. (2001), 2 
= ICorcoran &: Hamaguchil ((20071 ). 3 = [Parkin et al.. (.200a ). 4 = 
IPittard CorcoranI (|2002l ). 

plore the role of instabilities in the WCR, a grid-based 
hydrodynamic code is advan tageous for this pu rpose in 
comparison to an SPH code (|Agertz et alJ[2QQ7l ). 

2.3. X-ray emission 

To calculate the X-ray emission from the simulation 
we use emissivities for optically thin gas in collisional 
ionization equilibrium obtained from look-up tables cal- 
culated from the MEKAL plasma code containing 200 
logarithmically spaced energy bins in the range 0.1-10 
keV, and 101 logarithmically spaced temperature bins 
in the range 10"^ - 10^ K. When calculating the emer- 
gent flux we use energy depende nt opacities c alculated 
with version C 08. 00 of Cloudy (jFerlandl [2OOOI see also 
iFerland et al.lll998l ). The advected scalar is used to sep- 
arate the X-ray emission c ontributions from each w ind. 

Based on the findings of iParkin fc PittardI (|2010f ). we 
note that numerical heat conduction from the hot, low 
density postshock companion star wind to the cooler, 
higher density postshock gas of the primary star's wind 
will occur in the simulations. This effect is purely numer- 
ical in origin, and one must be careful when interpreting 
the derived X-ray emission. Therefore, as a cautionary 
measure we only consider the X-ray emission contributed 
by the postshock companion star's wind. 

3. THEORETICAL BACKGROUND OF RADIATIVELY 
DRIVEN WINDS 

3.1. Radiatively driven winds of single hot stars 

For the stellar atmosphere to be driven by the radia- 
tion field there is the simple requirement that GM/r^ < 
grad- The key to accurately describing the motion of the 
wind comes in the form of the ^rad term, and signifi- 
cant progress has been made using t he for ce multiplier, 
9Jt(t) = originally proposed bv IC AKI . 



where F is the radiative flux, and t = <Je'^thp[n • V(n • 
v)]~'^ is the Sobolev optical depth parameter. The force 
multiplier parameters a and k are determined from a 
power law fit to dJl{t). k can be interpreted as the fraction 
of flux that would be blocked if all lines were optically 
thick, a is the ratio of the line acceleration from optically 
thick lines to the total one. 

3.2. Finite disk correction factor 

A simplifying assumptio n adopted in early works on ra - 
diatively driven winds (e.g. lLucv fc SolomonlllQTOUCAKi ) 
was that of radially streaming photons, i.e. a point 
source of radiation. In this case the projected velocity 
gradient becomes dv^/dr where is the radial velocity, 
and angle integrals simplify to purely radial terms. Un- 
fortunately, the assumption of radially streaming pho- 
tons is a poor one close to the star where the wind is 
rapidly accelerated. To circumvent this problem it has 
become common place to incorporate a multiplicative 
factor to correct for the finite size of the stellar disk, 
which is attained by adopting the exact optical depth 
rather than the radial one. 



i^(r, dv/dr) \ 



grad (finite disk) 
grad (point source) 

(l+a)^+--(l + a/i;)^+- 
a(l+a)(l + a)-(l-/i2) 



(8) 



(9) 



.Fkt- 



grad 



where 6 is the angle subtended between the radial di- 
rection and a point on the stellar disk, jn = cos^, jn^ = 
VT^^C^^TO^, and a = {d\nv/d\nr) - 1 CCastod llQTOD . 
K is commonly ref erred to as the finite disk correc- 
tion f actor (FDCF) (jCastorlll974l : [CAH : [Pi[^ndrach et al.l 
Il986f ). 9Jt(t)FDCF is related to the standard (under the 
assumption of radially streaming photons) force multi- 
plier, 9Jt(t), via, 

3^(t)FDCF = Km{t) (10) 

In certain instances (e.g. iPauldrach et al.l 119861 : 
iStevens fc Pollo ck 1994) the FDCF is approximated as a 
purely radial, monotonic function (i.e. K{r^v^ dv/dr) ^ 
K(r)). Fig. [2] shows a comparison of monatonic and dy- 
namically consistent FDCFs; a monotonic FDCF always 
stays below 1, whereas a dynamically consistent FDCF 
rises above 1 in the wind acceleration region. For com- 
parison we have also plotted the numerical evaluation of 
the FDCF which we have implemented into FLASH - 
there is a good agreement with the analytical result. 

3.3. Interacting radiation fields in binary systems 

In the standard ICAKI model for radiatively driven 
winds the outwards acceleration is produced by the ra- 
diation field of the star. If one introduces an opposing 
radiation field two effects become apparent: radiative in- 
hibition ([ St evens fc Polloc k 19 9^) and radiative brakin g 
(jCavlev et al.l 119971 - see also lOwocki fc GavlevI [19951 ). 
The former refers to the reduction of the net rate of ac- 
celeration of a stellar wind due to the opposing radiation 
field, whereas the latter refers to the (sudden) decelera- 
tion of a stellar wind bef ore it reaches the WCR. In the 
IStevens fc Polloc3 (|1994f ) formulation of radiative inhi- 
(7) bition, the basic premise is that the line force (Eq. [7]) 
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100 150 200 250 

Distance from star (Rsol) 

Fig. 2. — Comparison of finite disc correction factors 
(FCDFs) calcul ated assuming monotonic flow (MONO, e.g. 
IStevens &l Pollocki i 19941, with dynamical consistence (DYN, i.e. 
including velocity gradi ent t erms - see Eq. [9]), and via a numerical 
evaluation (NUM - see ^ I2.1|) . The FDCFs are calculated using the 
companion (single) star wind velocity proflle (see Fig. O. 



becomes 



^rad 



{F1K1-F2K2} 



(11) 



where Fi, F2 are the radiative fluxes, and i^i, K2 are 
the finite disk correction factors for stars 1 and 2, respec- 
tively. The radiative inhib ition calculations presented by 
IStevens Pollockl (I1994D adopt ed a monotonic FDCF 
and, as noted by iGavlev et al.l (|1997[ ). it was for this 
reason that they did not see radiative braking in their 
calculations. 

Radiative inhibition is expected to be important if the 
stars are close enough that their opposing radiation fields 
reduce the net force within the acceleration region of 
the winds, whereas radiative braking is important in sys- 
tems where the imbalance in wind momenta is such that 
the stronger wind approaches the companion star close 
enough for its radiation pressure to exceed the wind ram 
pressure. Depending on the stellar and system param- 
eters adopted both effects may be important for 77 Car 
(|Parkin et al.l 120091 ). 

3.4. Considerations for r] Car 

The binary orbit of 77 Car is highly eccentric. There- 
fore, for a large part of the orbit the stars are well sep- 
arated, the winds have accelerated to terminal velocity 
when they reach the WCR, and the effects of radiative 
inhibition and radiative braking are negligible. However, 
an accurate description of wind acceleration (particularly 
of the companion's wind) is essential around periastron 
when the WCR enters into the wind acceleration regions. 

Compared to a monotonic FDCF, a dynamically con- 
sistent FDCF will provide additional acceleration close 
to the star driving the wind (see § 13. 2p . Therefore, one 
may wonder what difference using dynamically consis- 
tent FDCFs in Eq. [11] has on the level of radiative inhi- 
bition? In Fig.[3]we show the results from static two-star 
radiative driving calculations for r] Car at periastron. Us- 
ing a monotonic FDCF and with ki = k2 and ai = 0L2 
(i.e. the coupling between the primary's radiation and 
the companion's wind uses k = 0.50 and a = 0.68) we 
see that the terminal velocity attained by the compan- 



ion star wind is dramatically reduced from 3000 km s ^ 
(without RI) to 1500 km s"^ (RI - mono FDCF). If 
we include a dynamically consistent FDCF the level of 
radiative inhibition is lessened and the maximum veloc- 
ity attained now increases to 1600 km s~^ (RI - dynm 
FDCF). Interestingly, the velocity profile now displays 
radiative braking. Comparing the radiative driving cal- 
culations against those from a numerical hydrodynamic 
simulation calculated using the code described in § 12.11 
(HS), we see good agreement, albeit with a minor differ- 
ence in the position of the braking radius (the radius at 
which rapid deceleration of the wind occurs). 
In both the radiative inhibition fStevens & Pollock 



1994 ) and braking (Owocki & Gayley 1995; Gayley e t al.l 
19971 ) theories the line force is assumed to take the form 
of Eq. [m and in so doing the force multiplier param- 
eters used to accelerate/decelerate the stellar wind are 
assumed to be equivalent (i.e. k\ = k2 and ai = 0L2 
within the wind of star 1). At the alternative extreme, 
the line force could be described a^ 

fed = ^{9JTi(t)Fii^i - ^2(t)F2K2} (12) 

where 9}Tw(^) = k^t~^^ and the subscript w denotes the 
wind which the line driving parameters are describing. 
Fig. [3] shows the result of a hydrodynamic calculation 
(with radiatively driven winds) where ^rad is equivalent 
to that of Eq. [12] (i.e. the coupling between the primary's 
radiation and the companion's wind uses k = 0.30 and 
OL — 0.52). The resulting level of radiative inhibition 
is lower and the companion's wind reaches a maximum 
velocity of 2200 km s"^ (HS - A^i = A;2, o^i = ^2). Fur- 
thermore, the braking radius has increased by a factor of 
2. The explanation for this can be found by consider- 
ing the line driving parameters of each individual wind. 
The significantly higher luminosity of the primary star 
means that to achieve its huge mass-loss rate we actu- 
ally require lower values of k and a than the companion 
star (with its lower luminosity) requires to drive a weaker 
wind. In essence, the coupling between the primary's ra- 
diation field and its wind is weaker than the respective 
case for the companion star. Therefore, if the primary's 
radiation field is more strongly coupled to the compan- 
ion's wind (i.e k\ = /c2 and ol\ = 0^2), the decelerating 
force, and thus the net level of radiative inhibition, will 
be greater. Clearly, this result has implications for the 
dynamics of the WCR, and the infiuence o f radiative in- 
hibit ion on the X-ray emission from 77 Car (jParkin et al.l 
I2OO9I ) . In the following work we conservatively adopt the 
assumption of weaker coupling (Eq. [T2|) . but note that 
it could be stronger in reality. If this were the case one 
would expect lower preshock velocities for the compan- 
ion's wind, and reduced X-ray emission, particularly at 
higher energies (e.g. E > 2 keV). 

4. PERIASTRON SIMULATIONS 

As the stars move towards periastron passage their sep- 
aration contracts, their orbital velocities increase, and 

^ It is currently unclear what the correct choice for the cou- 
pling between the radiation field of a star with the opposing star's 
wind is. We merely note that previous works on the interaction 
of radiation fields i n binary sta r systems (e.g. Pittard et al. 1998; 
IGavlev et aP IMf: St-Lou is et al. 2005 ; Parkin et al. 2009) sug- 
gest that the actu al c oup ling lies somewhere between the two cases 
presented in Eas. Illlandll2l 
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Without Rl 
Rl - mono FDCF 
Rl - dynm FDCF 
HS 

HS - kH=kp, aH=ap 



50 100 150 200 250 300 350 

Distance from companion star (Rsol) 

Fig. 3. — Companion star wind velocity along the line of centres 
between the stars. The primary star is situated at a distance of 
359 R0 (corresponding to periastron separation when e = 0.9). The 
following solutions are shown: without radiative inhibition (RI)(i.e. 
single star), Rl using a monotic FDCF, Rl using a dynamically 
consistent FD CF, an d a numerical hydrodynamic simulation (HS) 
with separate ICAKI parameters for each wind, and a HS but with 
ki and ai used to drive both winds. The stellar and line driving 
parameters used to perform these calculations are noted in Tabled 



the preshock and postshock gas density (and velocity 
in the accelerating winds case) undergo a considerable 
rate of change as the WCR moves closer to the stars. 
To form a basis for qualitative and quantitative pre- 
dictions for the large-scale simulations that will follow 
(§ El), we begin with simulations focused on periastron 
passage where orbital motion is not included. A fixed 
resolution grid with cell size of c:^ 9.8 x 10^^ cm is used 
(i.e. no AMR) and the simulation box extends from 
X = z = ±2.34 X 10^^ cm and = - 1.88 x 10^^ cm 
(corresponding to xxyxz = 32 x 192 x 32 cells). The 
companion star is situated at (x, z) = (0,0,0) and the 
primary star is situated at (x, ^, z) = (0, 2.5 x 10^^ cm, 0). 
Zero-gradient boundary conditions are used on all faces 
of the simulation box. The calculations were evolved 
for a time of 4 x 10^ s, corresponding to 0.02 in or- 
bital phase, and therefore sufficiently long to examine the 
growth of instabilities over the characteristic time that 
the stars spend at closest approach. For comparison we 
have performed simulations with instantaneously acceler- 
ated winds (model Peri-IA) and radiatively driven winds 
(model Peri-RD). We now proceed with a discussion of 
the results of these models, including calculations of the 
intrinsic 7-10 keV X-ray luminosity and its variability, 
and also the results of additional tests performed with 
varying grid resolutions and at separations correspond- 
ing to orbital phases around periastron. 

The periastron simulations are initialized with two op- 
posing ffows separated by a contact discontinuity at the 
ram pressure balance point between the winds. As model 
Peri-IA is evolved through time the winds collide and 
shocks are formed. The primary's relatively slow and 
dense wind cools rapidly to form a thin dense layer, 
whereas the faster, more tenuous w ind of the compan- 
ion star behaves quasi-adi a baticallv (jPittard et al]|1998l : 
iPittard fc CorcoranI [2QQ1 iParkin et al.l l2009[ ). Conse- 
quently, there is a large density and temperature con- 
trast between gas on either side of the contact disconti- 
nuity and the growth of Ray leigh- Taylor (RT) instabili- 
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Fig. 4. — Density (upper row) and temperature (lower row) snap- 
shots of the X — y plane from models Peri-IA (left column) and 
Peri-RD (right column) at a simulation time of 2.4 x 10^ s. The 
companion star is situated at the lower boundary and the primary 
star is situated beyond the opposite end of the box. All plots show 
a region oi x = ±2.34 x lO^^ cm and y = (0 - 1.88 x lO^^ )cm - tick 
marks correspond to a distance of 1 x lO-*^^ cm. The hot postshock 
companion wind in model Peri-IA acts as a "cushion" against thin- 
shell instability. In contrast, NTSIs are permitted to grow by the 
lower temperature postshock companion's wind in model Peri-RD. 
A volume rendering of the model Peri-RD density distribution is 
shown in Fig. [S] 
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Fig. 5. — Density volume rendering of model Peri-RD at a sim- 
ulation time of 2.4 x 10^ s. The companion star is situated to the 
right of the image (the end face of the box cuts through its cen- 
tre) and the primary star is situated beyond the opposite end of 
the box such that only a small area of its surface fits within the 
box. Driven by instabilities in the WCR, fragments of cold, dense 
postshock primary wind gas (blue) can be seen colliding against 
the companion star (purple). 
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Fig. 6. — Variability of the intrinsic 7-10 keV X-ray luminosity 
for models Peri-IA and Peri-RD. Note that the simulation grid 
for models Peri-IA and Peri-RD is focused on the region between 
the stars, therefore the total X-ray emission from the entire WCR 
cannot be calculated. However, our comparison here highlights 
important differences between the flux levels of the two models. 



ties ensues. Subsequent perturbations to the co ntact dis- 
conti nuity seed the hnear thin-shell instability (jVishniad 
[1983"). However, the hot, quasi-adiabatic postshock com- 
panion's wind has a sufficiently fast sound speed, and 
thus rate of response to perturbations in the contact dis- 
continuity by the dense thin layer, to render it linearly 
stable and the re gime of non-lin ear instability growth 
is never reached (jVishniad 1 1994f ). Therefore, the ther- 
mal pressure of the hot postshock companion's wind 
acts like a cushion and prevents an instability-driven col- 
lapse/collision of the WCR against the companion star 
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Fig. 7. — Preshock companion wind speed measured along the 
line of centres as a function of orbital phase for simulations Peri-RD 
and Orbit-RD. 



(Fig.Sl). 

However, the inclusion of the radiative driving of the 
stellar winds in model Peri-RD, and thus of the wind 
acceleration regions, introduces some important differ- 
ences to this picture. Firstly, by a simulation time of 
t = 10^ s the primary's radiation field has reduced the 
preshock companion's wind velocity from 2500 km s~^ 
to 2200 km s~^, and radiative inhibition of the initial 
(single star) velocity profile (see Fig. [3j) is evident. As 
such, the postshock gas temperature is lower. Secondly, 
the WCR moves slightly closer to the companion star to 
attain a new ram pressure balance and in so doing the 
preshock (and thus postshock) gas densities are higher. 
The combination of these factors causes radiative cooling 
to become important for the postshock companion wind, 
and as the gas cools it loses the thermal pressure sup- 
port which ensures its stability. When this occurs NTSIs 
begin to grow and by a simulation time of 7 x 10^ s oscil- 
lations in the WCR push a dense clump of postshock pri- 
mary wind towards the companion star. This clump then 
fragments from the WCR and passes deep into the com- 
panion's wind acceleration region, narrowly missing the 
star. Driven by the NTSI, the position of the WCR con- 
tinues to oscillate throughout the simulation causing the 
mean preshock velocity of the companion's wind between 
the stars to vary in the range ~ 2100 — 2350 km s~^. At 
a time of 2.4 x 10^ s multiple fragmented clumps of post- 
shock primary wind collide against the companion star 
(see Figs. |4] and [5]). Such collisions occur at multiple 
times throughout the simulation, interspersed by brief 
recoveries as the shock oscillates away from the compan- 
ion star. 

An obvious question is why do collisions occur in model 
Peri-RD but not in model Peri-IA? As stated above, 
when radiative cooling becomes important for the post- 
shock companion's wind in model Peri-RD linear thin- 
shell instabilities can grow non-linearly with devastating 
results for the WCR. Hence, we can quantify the dif- 
ference between m odels Peri-RD and Peri-IA using the 
cooling parameter () Stevens et al.|[T992[ ) 



X 



^cool 
^flow 



^12 



M_ 



(13) 



where tcooi is the cooling time, tnow = 10"'^^(ii2/sps is the 
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flow time, Sps is the postshock sound speed along the line 
of centres, vs is the preshock gas velocity (in 10^ cm s~^), 
di2 is a characteristic distance (in 10^^ cm), and M_7 
is the mass-loss rate (in lO~^M0yr~^). Here we take 
di2 to be the distance of the companion shock from the 
star as we are interested in the importance of cooling 
close to the WCR apex. As a rule of thumb x > 1 
indicates quasi-adiabatic gas whereas X ^ 1 indicates 
that gas cools rapidly. Examining x foi" the compan- 
ion's wind we find values of 3.6 and 1 for models Peri-IA 
and Peri-RD, respectively. Therefore, on the timescales 
that we are interested in (i.e. the closest approach of 
the stars around periastron passage) the differences in 
postshock gas conditions between models Peri-RD and 
Peri-IA are sufficient to trigger rapid radiative cooling 
and the subsequent growth of catastrophic instabilities in 
the former. Interestingly, the differences between post- 
shock gas conditions in models Peri-RD and Peri-IA are 
relatively minor, but because of the strong dependence 
on the preshock velocity there is a fine line between a 
WCR which is only plagued by instabilities and one in 
which the instabilities ultimately destroy it. The differ- 
ences in the X-ray minima (extended in 1998 and 2003.5 
versus the quick recovery seen in 2009), may therefore re- 
flect differences in the behaviour of the WCR caused by 
a small change in one or both winds. We note again that 
our adopted radiation- wind coupling essentially provides 
a lower limit to the degree of radiative inhibition of the 
companion's wind. If a stronger coupling were adopted 
the companion wind would collide at lower speeds and 
radiate more efficiently, and a WCR disruption (poten- 
tially even more ferocious) as in model Peri-RD would 

o ccur more readily . 

lAkashi fc Sokerl ()2010l ) have recently presented numer- 
ical simulations akin to model Peri-IA. However, while 
the radiation fields were ignored, the gravitational influ- 
ence from only the companion star was included. This 
results in the balance point of the WCR shifting closer to 
the companion star, thereby increasing the cooling rate 
of the shocked gas from the companion's wind, and mak- 
ing it easier for instabilities in the WCR to strike the 
companion star. Clearly, one must include all relevant 
forces acting on the shocked gas. 

4.1. X-ray emission 

Previous models of energy dependent X-ray emission 
from T] Car have revealed that if the preshock compan- 
ion's wind is assumed to be at terminal velocity at or- 
bital phases coinciding with the extended X-ray min- 
imum there is an overestimate of roughly an order of 
magnitude in the 7-10 keV flux when com pared to that 
observed by XMM-Newton (|Parkin et al.l 12009). To ex- 
amine whether the differences between models Peri-IA 
and Peri-RD can provide clues to alleviating this dis- 
crepency we have performed X-ray calculations on the 
simulation output (Fig. [6]). The X-ray luminosity rapidly 
rises as the winds start to collide in the simulation and 
then reaches a peak before declining. In model Peri-RD 
the flrst minimum corresponds to the flrst occurence of 
a fragment of dense postshock primary star wind coming 
into close proximity and/or colliding with the compan- 
ion star. In contrast, no such collisions occur in model 
Peri-IA, and the variability in the 7-10 keV flux instead 



corresponds to oscillations in the shock front. The lower 
preshock velocities, and thus lower postshock gas temper- 
atures (see Fig. HI), in model Peri-RD combined with the 
shock obliquity created by the NTSI results in an average 
intrinsic 7-10 keV luminosity which is roughly a factor of 
8 lower than that calculated for model Peri-IA. There- 
fore, it seems that a substantial disruption of the WCR 
can provide a plausible explanation for the observed be- 
haviour of the hard X-ray flux (Hamaguchi et al..,20Q7l ). 

4.2. Resolution dependence 

The growth rate of instabilities are dependent on the 
wavelength, with smaller wavelengths producing faster 
growth rates. To examine the dependence on grid 
resolution we have performed testfl with cell sizes of 
(4.9 — 39.2) X 10^^ cm. Encouragingly, the general charac- 
teristics appear consistent. In each case instabilities grow 
and break-up the WCR causing stochastic collisions be- 
tween dense clumps and the companion star. There does 
not appear to be any correlation between the grid reso- 
lution and the frequency of collisions. 

Additional tests were performed using simulations 
boxes with a cell size of 9.8 x 10^^ cm and a width up 
to X = z = ±3.5 X 10^^ cm. Consistent with model Peri- 
RD, when the acceleration (and radiative inhibition) of 
the companion's wind are considered the WCR becomes 
disrupted. Furthermore, the signiflcantly lower 7-10 keV 
X-ray flux illustrated by the comparison between models 
Peri-RD and Peri-IA is reproduced. 

4.3. Variation of preshock velocity with orbital phase 

The separation of the stars changes rapidly around pe- 
riastron passage and, therefore, so does the level of radia- 
tive inhibition, and the preshock companion star wind ve- 
locity. To examine this orbital phase dependence we have 
repeated model Peri-RD with stellar separations corre- 
sponding to orbital phases in the range (j) = 0.96 — 1.04, 
the results of which are plotted in Fig. [71 Error bars 
have been included to indicate the range of preshock ve- 
locities attained due to oscillations in the position of the 
companion shock. As the stellar separation tends to- 
wards its periastron value the instabilities become more 
vigorous, perturbing the shock front in the wind accel- 
eration region (where the velocity gradient is high) and 
thus broadening the range of sampled velo cities. 

Com paring Fig. [7l against flgure 21 of iParkin et al.l 
(|2009l ) it is evident that the conservative assumption of 
weak coupling adapted in the present paper leads to a re- 
duced level of radiative inhibition and noticeably higher 
preshock velocities (even though fragments of the WCR 
stochastically collide with the companion star at (/) = 1.0 
an average preshock velocity of 2200 km s~^ is at- 
tained jj. 

^ Note that to accurately model the acceleration of the stellar 
wind requires that the resolution not decrease in regions of sharp 
velocity gradient (i.e. within a radius of ~ 3 R^). Hence in the 
resolution tests the companion star wind acceleration region is cov- 
ered with cells of side ^ 9.8 x 10-'^'-' cm and we only vary the cell 
size for postshock gas using the AMR capability of the code. 

^ In this work the winds are initiated by forcing an outflow in a 
spherical shell around each star. This prevents accretion from over- 
powering the wind, which may s hut down the companion's wind 
for a signiflcant pe riod of time (|Sokeii 120051 : lAkashi et al.l 120061 : 
IKashi So'ra[2Q09al ) 
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5. ORBIT SIMULATIONS 

While the periastron simulations presented in § [H pro- 
vide some interesting insight, to fully understand the flow 
dynamics and resulting X-ray emission from rj Car we 
must perform simulations which include the motion of 
the stars. For this purpose we have performed two large- 
scale simulations: one in which the winds are assumed to 
be instantaneously accelerated at the surface of the star 
(model Orbit-IA), and another where the winds are ra- 
diatively driven (model Orbit-RD). Our adopted orbital 
and stellar parameters are noted in Tables [1] and [2l The 
mapping of the winds in model Orbit-IA is described in 
Appendix \^ 

To model the orbit of the r] Car binary system we 
use a simulation domain which extends from x = y = 
±2 X 10^^ cm and z = (0 - 2 X 10^^)cm. Outflow bound- 
ary conditions are used for all faces of the simulation box 
except the lower z boundary which is a symmetry bound- 
ary. The grid is initialized with xxyxz = 16x16x8 cubic 
blocks each containing 8^ cells. We allow for 9 levels of 
refinement, which results in an effective resolution on the 
finest grid level ofxxyxz = 32768 x 32768 x 16384 (i.e. 
a cell size of 1.22 x 10^^ cm or 1.75 R©). The refinement 
of the ^rid is deterni ined using a second-derivation error 
check (jFrvxell et al.| [2QQQ) on p and the requirement of 
an effective number of cells between the stars to accu- 
rately describe the WCR dynamics (see Appendix [B]) . 
The former identifies cells for refinement (and derefine- 
ment), whereas the latter controls the maximum resolu- 
tion (for postshock gas) at a given orbital phase. 

In the following we discuss the dynamics from sim- 
ulations Orbit-IA and Orbit-RD, focusing first on the 
large-scale dynamics and then concentrating on the oc- 
currences around periastron passage. We then present 
the results of X-ray calculations performed on the simu- 
lation output and detailed comparisons against observed 
lightcurves and spectra obtained with RXTE and XMM- 
Newton. 

5.1. Dynamics 

5.1.1. Large-scale dynamics 

For most of the orbit the stars advance relatively slowly 
and the effects of orbital motion on the gas dynamics 
are modest. At these times the separation of the stars 
is sufficiently large for the stellar winds to be at ter- 
minal velocity when they reach their respective shock. 
However, around periastron passage the orbital veloci- 
ties and stellar separation change rapidly as the stars 
career past each other. The rapid motion of the stars 
around periastron passage acts to contort the WCR into 
a spiral-like sh ape which subsequently expands ou twards 
with the flow (jOkazaki et al.l I2QQ8I : iParkin et aP [20091. 
Figs. [8] and [9] show a series of gas density and temper- 
ature snapshots from simulations Orbit-IA and Orbit- 
RD, respectively, from which the salient features of the 
wind-wind collision and the effects of orbital motion can 
be seen. The contact discontinuity separates the winds, 
being abutted by a thin dense layer of postshock pri- 
mary wind on one side, and a thicker "puffed- up" region 
of high temperature quasi-adiabatic postshock compan- 
ion wind on the other. The density contrast and veloc- 
ity shear across the contact discontinuity subject it to 
Kelvin-Helmholtz (KH) and RT instabilities. 



At (j) = 0.5, the stars are at their largest sepa- 
ration, moving at their slowest orbital velocities, and 
it is apparent that a reasonable approximation to 
the X-ray emitting part of the WCR can be at- 
tained from a 2D model w hich neglects orbital motion 
(jPittard fc CorcoranI 120021 ) . The apex of the WCR is 
closer to the primary star in model Orbit-RD than in 
model Orbit-IA due to the primary's wind ram pressure 
being slightly lower at an equivalent radius in the former 
- this is a numerical artifact resulting from the differ- 
ent ways by which the winds are initiated in the models. 
By (j) = 0.9 the stars have moved around in their orbits 
and their separation has reduced, yet the general shape 
of the large scale WCR has not changed considerably. 
However, the rapid motion of the stars during a brief 
period around periastron causes the WCR close to the 
stars to become very distorted (see Figs. [T2l and [T3|). al- 
though the large scale WCR is not impacted upon by 
periastron passage at = 1.0. At (/) = 1.1 the contortion 
of the WCR has had time to advect out in the flow, and 
the spiral structure of the WCR visible in previous 3D 
mode ls is reproduced (jOkazaki et al.|[2QQ8l : iParkin et al.l 
l2009f ). Interestingly, the increasing influence of orbital 
motion as the stars approach periastron results in an 
increasingly asymmetric temperature distribution in the 
postshock gas, with gas in the leading arm of the WCR 
clearly at a higher temperature than gas in the trailing 
arm (Figs. [8]and[9j). In model Orbit-IA the postshock 
gas in the trailing arm is far from smooth; there are thin, 
dense fllaments separated by low density, rarefled gas. 

Fig. Ho] displays a snapshot of the 3D structure of 
the WCR at periastron in simulation Orbit-RD. The 
clumpy /rippled texture to the WCR is clearly evident, 
as is the immense difference in scale between the large 
scale structure of the WCR and the region between the 
stars where the most ferocious activity occurs. 

The snapshots of simulations Orbit-IA and Orbit-RD 
in Figs. [8] and [9] re veal a number of noticeable differences. 
Firstly, the opening angles of the primary star and com- 
panion star shocks are slightly larger in model Orbit- 
RD when compared to model Orbit-IA. This occurs due 
to the radiation flelds of the stars; off- axis the preshock 
speed is greater than the on-axis preshock speed. The 
radiation fields of the stars also cause differences in the 
postshock gas between the two simulations. For instance, 
there are less instabilities in the WCR in model Orbit-RD 
when compared to model Orbit-IA. Examining the layer 
of postshock primary wind in the two models we find 
that the gas velocity (density) is higher (lower) in model 
Orbit-RD. Noting that the growth rate for KH instability 
is oc a/piP2/{pi + P2), the lower density and higher ve- 
locity in model Orbit-RD results in slower growth rates. 
Therefore, the fiow exits the grid with fewer noticeable 
perturbations. But why does this difference occur? Con- 
sider the motion of postshock gas away from the stag- 
nation point. In model Orbit-IA postshock gas is ac- 
celerated by the pressure gradient. However, in model 
Orbit-RD the postshock primary wind is of a low enough 
temperature to be acted upon by the radiation fields of 
the stars. Therefore, the higher velocity of postshock 
primary's wind in model Orbit-RD is due to radiative 
acceleration. 

Examining the leading arm of the WCR the layer of 
postshock primary wind appears thicker in model Orbit- 
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Fig. 8. — Snapshots of the gas density (left column) and tem- 
perature (right column) in the orbital {x — y) plane from model 
Orbit-IA at (/) = 0.5 (top row), 0.9 (upper middle row), 1.0 (lower 
middle row), and 1.1 (bottom row). The orbital motion of the stars 
is calculated in the centre of mass frame. At apastron (0 = 0.5) 
the primary star is to the right, and the companion star is to the 
left, of the image centre. The motion of the stars proceeds in an 
anti-clockwise direction. All plots show a region of ±2 x 10-*^^ cm - 
large axis tick marks correspond to a distance of 1 x 10^^ cm. 



RD when compared to model Orbit-IA. Interestingly, at 
phases close to periastron when this gas resides close to 
the stars, its thermal pressure is lower than the radia- 
tion pressure, which provides resistance against contrac- 
tion and thus widens the layer. However, comparing the 
snapshots at = 1.1 we see that at an equivalent dis- 
tance from the stars the density of the unshocked winds 
is slightly higher in model Orbit-RD, which means the 
mass in the swept-up shell is greater. The inertia of the 
swept-up mass is therefore greater in model Orbit-RD, 




Fig. 9. — Same as Fig. [8] except model Orbit-RD is shown. 



which accounts for the smaller distance that the spiral 
has travelled to by (/) = 1.1 in this model compared to 
model Orbit-IA. 

The width of the dense layer clearly affects the growth 
of instabilities in the expanding spiral-shaped shell - in 
model Orbit-IA the shell appears to be subject to the 
NTSI, whereas in model Orbit-RD the additional thick- 
ness to the layer renders it stable. This is unsurpris- 
ing as the stability of an expanding shell depends o n 
the shock thickness (jVishniadllQSllWiinsch et al.ll2QlQ[ ). 
This raises questions about the fate of the expanding 
shell in each simulation. As its outwards acceleration 
is decreased by an increasing amount of swept up mass 
its Mach number will decrease and as the shocks dissi- 
pate it will gradually mix with the bubble of companion 
wind which it encases. However, this only appears to be 
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Fig. 10. — 3D gas density volume rendering of simulation Orbit-RD at periastron passage (cf) = 1.0) viewed from above (left panel) and 
below (right panel) the simulation box. The entire simulation domain is shown. 




Fig. 11. — Snapshots of orbital {x — y) plane from simulation 
Orbit-RD (left panel) and Orbit-IA (right panel) at (/> = 1.1 show- 
ing the fluid dye variable. Values of 1 and correspond to a cell 
consisting entirely of primary or companion wind material, respec- 
tively, and intermediate values indicate a mixture. The plots show 
a region of ±2 x 10-*^^ cm - large axis tick marks correspond to a 
distance of 1 x 10-*^^ cm. For the corresponding density and tem- 
perature images see Figs. [8] and [O] 



the case for model Orbit-RD in certain directions. In the 
trailing arm of the WCR the dense layer of postshock pri- 
mary wind is photo- ablated by the radiation fields. Yet, 
this only occurs in the trailing arm of the WCR because 
the dense layer can be driven into a more tenuous cav- 
ity made by the companion's wind as it swings through 
periastron. In contrast, the leading arm is bordered by 
unshocked primary wind, and therefore the dense layer 
is merely widened by the radiation pressure. Examining 
the fluid dye variable (which tracks the quantity of each 
wind in a given cell) the photo-ablation by the radiation 
fields of the stars appears to be a very effective mixing 
agent and (unlike in model Orbit-IA) smooth out any 
filamentary structure. Comparing models Orbit-RD and 
Orbit-IA at = 1.1 we see that in the former the photo- 
ablation of the dense layer snips the tail of companion 
wind gas in the WCR (Fig. [11]). 

Following the expansion of the shell as it exits the sys- 
tem, we see that the shocked gas in the current spiral 
crashes through the remnant of the WCR from prior 
to periastron passage, and Richtmyer-Meshkov (RM) in- 



stabilities are formed at the interface between the two. 
The additional mixing provided by photo-ablation and 
the collision of successive expanding spiral-like shells will 
have implications for models of observations directly re- 
lated to the posts hock primary st ar's wind, e.g. forbid- 
den l ine emission (jGull et al.ll2QQ9f ). episodic dust forma- 
tion (lSmithll2QlQ[). and high- velocity absorption features 
(|Groh et al.ll2010b[ ). 

5.1.2. Periastron passage 

Focusing now on the periastron passage of r] Car, where 
there is a considerable rate of change in the position of 
the stars and the distribution of the winds, a number of 
interesting features can be identified. Firstly, there ap- 
pears to be a sudden change in the level of instability 
in the shocks as periastron is approached. For instance, 
prior to periastron the shocks seem to be unstable in both 
model Orbit-IA and Orbit-RD (see cj) = 0.990, and 0.995 
snapshots in Figs. [T2l and [T3|) . and one sees that, consis- 
tent with the predictions of the periastron simulations in 
§ m the postshock gas is more unstable when the wind 
acceleration regions are taken into consideration (model 
Orbit-RD) than in the case where they are not (model 
Orbit-IA). This fact is highlighted at = 0.995 where 
in model Orbit-RD a shard of dense postshock primary 
wind gas is forced through the shocks by instabilities in 
the leading arm of the WCR and comes close to colliding 
against the companion star before being rapidly ablated 
by its wind. However, as = 1.000 is approached, and 
the orbital velocities of the stars increase, the shocks di- 
rectly between the stars appear much more stable and the 
catastrophic disruption of the WCR seen in model Peri- 
RD is not reproduced in model Orbit-RD. The reason 
for this is two-fold. Firstly, the level of radiative inhibi- 
tion of the preshock companion's wind is not as severe 
when orbital motion is included - this point is discussed 
in more detail below. And secondly, orbital motion, and 
the velocity component that it provides perpendicular to 
the line of centres, increases the preshock velocity. Con- 
sequently, the postshock gas is less radiative and the pre- 
vailing thermal pressure stabilizes the WCR against the 
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Fig. 12. — Snapshots of the orbital (x - y) plane from model Orbit-IA showing (from top to bottom) </> = 0.990, 0.995, 1.000, 1.010, 
1.020, and 1.040, and (from left to right) density, temperature, and speed. At periastron (0 = 1.000) the primary star is to the left, and 
the companion star is to the right, of the image centre. The plots at </) = 0.990, 0.995, 1.000, and 1.010 show a region of ±1 x 10-*^^ cm, 
whereas the plots at = 1.020 and 1.040 show a region of ±5 x 10-*^^ cm. In all plots large axis tick marks correspond to a distance of 
±1 X 10^^ cm. Note the difference in colour scale used for the density plots at 6 =1.020 and 1.040. 
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Fig. 13. — Same as Fig. 1121 except model Orbit-RD is shown. 
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growth of NTSIs. These results imply that a significant 
disruption to the WCR must be afoot prior to the rapid 
rotation of the stars (i.e. prior to (j) = 1.0). 

At orbital phases after periastron the contortion of the 
WCR by the motion of the stars is most evident. The 
snapshots at = 1.01 and 1.02 in Figs. [T2l and [T3l show 
how the unshocked companion wind is completely en- 
closed by the WCR. Following periastron one can see that 
the arms of the WCR become twisted to such an extent 
that they collide downstream (see = 1.02 in Figs. [T2l 
and [13]), so that postshock gas in the trailing arm is fur- 
ther heated as it passes through the shocks in the leading 
arm. The wind of the companion star (on the side facing 
away from the primary star) now collides against the in- 
side of the leading arm of the WCR and helps to heat this 
gas. In model Orbit-RD the preshock companion wind 
in this direction is accelerated by the cumulative radia- 
tive driving force from both stars, the result of which is 
that immediately preshock this gas has reached a veloc- 
ity of 3300 km s~^, noteably higher than the (terminal) 
velocity of the companion wind in model Orbit-IA. Con- 
sequently, the ram pressure of the preshock gas is higher 
and it attains a pressure balance with the postshock gas 
at a larger distance from the companion star, causing 
the region of postshock companion wind to be thinner 
in model Orbit-RD compared to model Orbit-IA. How- 
ever, despite the higher ram pressure of the companion's 
wind in this direction, the expansion of the spiral at this 
time occurs at roughly the same speed because the com- 
panion's wind is eflFectively running into a wall of dense 
postshock primary wind with a high inertia. Therefore, 
it is the postshock primary wind which largely dictates 
the rate of expansion. Examining the tail of the WCR in 
the density snapshots at ^ = 1.01, 1.02, and 1.04 one can 
see a striking example of photo-ablation by the radiation 
fields of the stars (see Fig. [T3]) . 

5.1.3. Preshock velocities 

The simulations presented in § HI in which the stars 
were static with a separation corresponding to perias- 
tron, predicted that the winds shock at such speeds that 
the WCR becomes very unstable, so that dense clumps 
of gas may collide with the companion star. Yet, when 
we add orbital motion to the mix we do not see a disrup- 
tion of such magnitude. If anything, the rapid motion 
of the stars around periastron passage appears to have a 
stabilizing effect on the apex of the WCR. We can begin 
to dissentangle this puzzle by examining the preshock 
velocities from models Orbit-RD and Peri-RD (Fig. [7]). 
Clearly, at orbital phases leading up to periastron the 
companion's wind is not reduced to the same preshock 
velocity in model Orbit-RD as in model Peri-RD. In fact, 
in model Orbit-RD the preshock companion's wind ve- 
locity does not decrease to that in model Peri-RD until 
(j) 1.005, at which point the lowest preshock velocity 
is attained. Following this phase the agreement between 
model Orbit-RD and Peri-RD improves and there ap- 
pears to be a good correlation between the level of ra- 
diative inhibition occuring. This behaviour is due to the 
combination of two effects introduced by orbital motion 
(see Fig. [T4|) . In the frame of the WCR, the preshock 
gas velocity is equal to the wind velocity (in the frame 
of the star) plus the relative velocity of the star along 
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Fig. 14. — The orbital speed of the companion star around pe- 
riastron passage. The components for the transverse and hne-of- 
centres velocity, as well as the sum (total) are shown. For compar- 
ison, we also plot the difference between the preshock companion 
star wind velocities from models Orbit-RD and Peri-RD. Note that 
these values are calculated in the frame of reference of the centre 
of mass. 



the line of centres, (isep- Therefore, prior io (j) c:^ 0.99, 
dsep is positive and the contraction of the stellar separa- 
tion will increase the preshock wind velocity. However, 
dsep is positive until ^ ~ 0.99 and then becomes negative 
(Fig. [T4|) . so at > 0.99 the radial motion will decrease 
the preshock wind speed. The wind material also has 
a considerable velocity in the transverse direction which 
increases (decreases) as the stars approach (depart from) 
periastron. This motion reduces the wind density along 
the line of centres, and as ^rad the wind accel- 

eration increases, and the impact of radiative inhibition 
decreases (see Fig. [15]). Note that the alteration to ^rad 
will apply for both radiation fields, but since 0^2 > ai the 
reduction in p benefits the secondary's radiative driving 
more than inhibition by the primary. 

Prior to (j) 1.01 the combination of these two effects 
increases the companion's preshock wind speed along 
the line of centres compared to the static-star case (see 
Figs. [3l and [T4|) . However, following periastron the en- 
hanced preshock velocity introduced by transverse mo- 
tion is countered by the receding relative motion of the 
companion star along the line of centres. Examining the 
sum of the orbital velocity components around periastron 
one sees good qualitative agreement with the difference 
between models Orbit-RD and Peri-RD (Fig. [M]). We 
note that although the critical point radius of the wind 
acceleration is not resolved in the simulations, the infiu- 
ence of orbital acceleration on the preshock winds will be 
negligible (Appendix [U]). 

In simulation Orbit-IA a maximum wind speed of 
~ 3200 km s~^ is attained. When the acceleration of the 
winds by the radiation fields of both stars is considered 
the maximum speed for gas grows to ~ 4500 km s~^, 
with the highest velocities possessed by unshocked com- 
panion wind gas which has been subjected to the cumula- 
tive driving force from both stars. Although velocities of 
this order will produce extremely hot gas when it shocks 
it is unlikely that there will be a significant observable 
signature as the total gas mass, its postshock density, 
and thus the free-free emission, is very low. 
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Fig. 15. — An illustrative example of the effect of orbital motion 
on the resultant line of centres wind acceleration. When the stars 
are static the line of centres (dot-dashed line) is aligned with radial 
flow (left). However, when a star undergoes orbital motion (right) 
the flow along the line of centres at time t2 consists of wind material 
emitted at earlier times, ti. Consequently, the path length for 
flow to reach the line of centres is larger (i.e. R2 > Ri) and as 
p oc and gfj-ad oc p~^, this leads to an increase in the radiative 
acceleration. 



5.1.4. Resolution dependence 

Resolution tests were performed for the large-scale sim- 
ulations using finest cell sizes (for the postshock gas) 
in the range 0.66 - 2.44 x 10^^ cm (0.88 - 3.5 R©). 
As the simulation resolution is lowered the growth of 
KH instabilities at the contact discontinuity is less pro- 
nounced, and vice-versa for higher resolution. Impor- 
tantly, there was a consistent finding that a substantial 
disruption/collapse of the WCR does not occur in simula- 
tions which include orbital motion (for our adopted stel- 
lar, system, and radiation- wind coupling parameters). 

We note again that the key reason that a catastrophic 
disruption/collapse of the WCR does not occur at perias- 
tron in model Orbit-RD is that the companion's preshock 
wind velocity does not reach sufficiently low values to 
trigger radiative cooling in the postshock gas. Therefore, 
although the simulations are only progressed through a 
single orbit, continuing the simulation for subsequent or- 
bits would not result in a disruption/collapse. 

5.2. X-ray emission 

T] Car is an exceptionally bright source of X-ray emis- 
sion. This X-ray emission provides a direct probe of the 
WCR, and detailed comparisons of model predictions 
against observations can determine key parameters 
of the stars and their winds. To date, considerable 
focus has been placed on constraining the orientation 
of the orbit on the plane of the sky. There remains 
a lack of consensus on this issue, with some models 
preferring the companion star to pass behind the 
prim a ry star at periastron (iDa mineli 1996; Pittard et al.l 
1998; Corcoran et al. 2001; Pittard & Corcoran 2002"; 
Corcoran 2005; Akashi et al. 2006: Hamaguchi et al.. 



20071: iNielsen et al.l |2QQ7|; HBenlevet al.l |2QQ8 



Qkazaki et all 120081: [Parkin et al] 120091: Ipuh et al.l 



2(^1 iRichardson etldl [2M0I : IGroh et al.l [2of?)b[ ) 

whereas others prefer the orientation the opposite way 
around (^Falceta-Goncalves et al. 2005; Abraham et aL, 
2QQ5a|; [Abraham Falceta-Goncalv esI 1200?!: 

Kashi fc Sokerl 120081 : iFalceta-Goncalves fc AbrahamI 



I2OO9I : lAbraham fc Falceta-Goncalves"20 1Qh. or the sys - 
tem at quadrature (Ishibashi 2001; Sm ith et al.ll2004 ). 
The inclination of the orbital plane, and whether or not 
it is aligned with the equatorial skirt of the Homun culus 
nebula, is still debated (see iParkin et al.l I2009L and 
references there-in). Our new models can also be used 
to test the orientation. We define our line of sight 
geometry as follows: the inclination angle, i, is measured 
against the 2;-axis {i = 0° would view the system from 
directly above the orbital plane), and the angle is 
measured against the negative x-axis (companion star in 
front at apastron) such that 6 increases in the prograde 
direction {0 = 90° would align the line of sight with the 
negative 7/-axis). We adopt viewing angles of i = 42° 
and = 20°, in agreement with values determined 
from the mos t recent modelling o f the X-ray emissio n 
from T] Car (jOkazaki et al.l 120081 : iParkin et all l2009f ). 
In the following analysis we define (j) = 1.0, 2.0, and 
3.0 as periastron passage in 1998, 2003.5, and 2009, 
respectively. 

5.2.1. X-ray lightcurves 

The X-ray lightcurves for models Orbit-IA and Orbit- 
RD are shown in Fig. [1 6} The characteristic shape of 
77 Car's X-ray li^htcurv e (ICorcoran et al.ll2001l : [Corcoran] 
2005; Corcoran et al. 201Q[) is reproduced by the models; 
the X-ray fiux remains relatively fiat for the majority of 
the orbit {(j) 0.05 — 0.7) and then as the stars approach 
periastron there is a sharp increase in X-ray fiux (0 
0.7 — 0.99) followed by an abrupt decline to a minimum at 
periastron, and a recovery of the X-ray fiux as the stars 
move away from periastron. The difference in apastron 
ffux between models Orbit-RD and Orbit-IA is due to 
the WCR being closer to the primary star in the former 
(see § 15. ip . Due to the slightly different wind densities 
at equivalent distances from the stars in models Orbit- 
IA and Orbit-RD we do not attempt to infer suggested 
alterations to model parameters which would improve 
the model fits. However, we note that on the whole the 
adopted wind momentum ratio provides good agreement 
between the model and the observed lightcurve. 

Examining the orbital phase range (f) = 0.9 — 1.1 reveals 
some noteable differences between the observations and 
the models. Both model Orbit-IA and, to a lesser ex- 
tent, model Orbit-RD overestimate the peak in average 
ffux observed by RXTE prior to periastron. At (j) = 1.0 
the models do not reach the same minimum in ffux as 
the RXTE data. This is, at ffrst, a somewhat puzzling 
result as previous models were able to reproduce the low 
ffux lev el at = 1.0 via an eclipse of the X-ray emission 
region (! Okazaki et all 120081 : IParkin et all [2009'). How- 
ever, inspection of the 2-5 keV broadband image at this 
phase reveals considerable emission from downstream gas 
in the tail of the WCR (Fig. [T7|). This emission predom- 
inantly comes from the trailing arm of the WCR with 
a weaker contribution from the downtream gas in the 
leading arm (see Fig. [18] - note that due to our adopted 
viewing angle the emission close to the apex in the lead- 
ing arm is strongly absorbed by the intervening primary 
wind). This lack of agreement provides an important 
clue as to the nature of the X-ray minimum of 77 Car - 
to match the low ffux level at (j) = 1.0 the apex of the 
WCR must be disrupted at a slightly earlier phase such 
that the disruption propagates out and removes the X- 
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ray emitting plasma in the trailing arm. The residual 
flux observed by XMM-Newtoi^ at = 1.0 could then 
be contributed by the X-ray emitting plasma in the lead- 
ing arm of the WCR (i.e. in the downst ream gas behind 
the companion star - iParkin et al.|[2QQ9[ ). The flow time 
for the postshock companion star wind to be accelerated 
from the shock apex to point "A" in Fig.[TH]is 12 days 
(S(j) :^ 0.006), which suggests that the disruption of the 
WCR must be afoot by (/) = 0.994, in agreement with the 
0-40 days proposed bv iSokerl (|2005D for the commence - 
ment of an accretion event (see also lAkashi et al.]l2006f ). 

A key failing is that neither model reproduces the ob- 
served extended X-ray minimum. The X-ray flux from 
model Orbit-IA recovers hy (j) = 1.017, whereas model 
Orbit-RD shows a slightly extended minimum with a 
recovery hy (j) = 1.025. It would therefore seem that 
the width of the extended X-ray minimum cannot be re- 
produced when the spatial extent of the emitting region 
is taken into consideration, which poses difficulties for 
models based o n an e clipse of the X-ray emission region 
([Okazaki et I2008L see also Parkin et a l. 2009). A 
step is seen in the minimum of model Orbit-RD, which 
corresponds to lines of sight passing through the dense 
shell of postshock primary wi nd, although it is n ot as 
pronounced as in the model of iParkin et all (j2009[ ). En- 
couragingly, the width of the X-ray minimum in model 
Orbit-RD agrees well with the less extended cycle 3 mini- 
mum observed by RXTE in 2009 (see lower right panel of 
Fig. [16]). This result implies that in the previous two cy- 
cles the extended X-ray minimum resulted from a catas- 
trophic disruption of the WCR, whereas the shorter du- 
ration of the recent minimum indicates that the WCR 
was not so severely disrupted the last time around. 

5.2.2. X-ray line profiles 

If the apex of the WCR is heavily disrupted the pro- 
jected velocity along the lines-of-sight from high temper- 
ature. X-ray emitting postshock gas will be different. A 
high spectral resolution observation with Chandra found 
that at = 1.009 there was a positive line shift (i.e. away 
from the observer) for the S XVI line, w hereas S XV, 
Si XIV, and Si XIII had negative line shifts (jHenlev et al.l 
[1)08, see also Behar et al. 2007). In the absence of a sta- 
ble WCR between the stars the higher excitation S XVI 
line could still have been emitted by hot postshock com- 
panion star wind close to the WCR disruption which (for 
our adopted viewing angle) is travelling away from the 
observer (labelled "B" in Fig. [p). The S XV, Si XIV, 
and Si XIII lines could then originate from lower temper- 
ature gas further downstream in the leading arm of the 
WCR (labelled "C" in Fig. [18]), which is moving towards 
the observer, albeit at a steep angle, such that there will 
be a negative projected hue shift. 

5.2.3. Column density 

The emission- weighted column densitie^ from models 
Orbit-IA and Orbit-RD appear to be largely similar for 

^ During the deep minimum there will be some contamination 
from the cosmic background in the RXTE data. 

^ The emission-weighted column density (EWC) is calculated as 
Newc = ^^H^intx/^^intx, where A^H and Lintx are the column 
density and 1-10 keV intrinsic luminosity from a given line of sight, 
and the summation is over all sight lines (pixels) in the X-ray 
image. 



the majority of the orbit (Fig. [19]), though model Orbit- 
IA shows a slower decline in column density in the orbital 
phase range 1.0-1.4. This difference is caused by the dis- 
tribution of the dense layer of postshock primary wind 
in relation to the lines of sight to the X-ray emitting 
plasma (see Figs. [H [El and [H]) . In model Orbit-IA 
the growth of KH instabilities at the contact discontinu- 
ity in the leading arm of the WCR distort the dense layer 
of postshock primary wind and push it into the path of 
X-rays as they exit the system, consequently increasing 
the column density. In contrast, the lack of KH instabil- 
ities in model Orbit-RD results in a more rapid decline 
in column density as the stars move away from perias- 
tron. Interestingly, the column density calculated from 
model Orbit-RD ag rees very wel l with the results from 
the model used bv IParkin et al.l ()2QQ9f ). which did not 
allow dynamical instabilities. 

The column densities derived from fits to XMM- 
Newton spectra in the 2-10 keV energy range 
(jHamaguchi et al"]l2QQ7[ ) are al so plotted in Fig. [191 At 
(j) 1.47 the models and the iHamaguchi et all ()2QQ7l ) 
values are in good agreement, with a rather modest re- 
duction of 1.15 required to align the models and obser- 
vation. Our adopted viewing angle places the compan- 
ion star in front of the primary at this orbital phase. 
Therefore, Nn oc M2, and the aforementioned reduc- 
tion factor would bring the companion star mass-loss 
rate to a value of 1.2 x lO~^M0yr~^. However, at 
phases close to periastron the column densities calcu- 
lated from models Orb it-IA and O rbit-RD are noticeably 
higher than the Hama guchi et al.l (|^07) values. During 
this part of the orbit there is a more complicated distri- 
bution of both the shocked and unshocked winds which 
lines of sight to the X-ray emitting plasma must pass 
through. Therefore, we cannot infer such a clear-cut re- 
duction in mass-lo ss rates from a compar ison between 
the models and the ^Ha maguchi et al.l (|2QQ7l ) column den- 
sities. At face value. Fig. [19] implies that a reduction in 
the mass-loss rates of at least a factor of 2 is necessary, 
which gives lower limits to the mass-loss rates of Mi = 
2.4 X 10-"^ M© yr-i and M2 = 7 x 10'^ M© yr-^ . How- 
ever, this seems unlikely, given other evide nce for mass- 
loss rates g r eater than those assum e d here (jHillier et al.l 
I2QQII . [20061 : Ivan Boekel eFall I200I iGroh etani2010a[ ). 
unless these estimates have in turn been "contaminated" 
by emission from the WCR. Another, perhaps more 
palatable way to reconcile the difference would be to in- 
voke the presence of (X-ray) ionized gas close to the stars, 
which would reduce the opacity and thus allow for higher 
mass- loss rates. Alternatively, perhaps observationally 
derived column densities around periastron are suspect, 
due to the complexities of disentangling different spec- 
tral components, and/or the vagaries o f fittings; simpler 
plasma models to complicate d systems (jAntokhin et al.l 
120041 : iPittard Parkinll2010D . For instance, any resid- 
ual soft X-ray emission near r] Car (unresolved to XMM- 
Newton) will reduce the apparent column density. 

Despite the model column densities being higher than 
the observed values, we can extract some useful infor- 
mation from a comparison of the variation in column 
density as a fun ction of orbital phase. F or instance, the 
models and the IHamaguchi et al.l (|2007[ ) values agree in 
that there is a sharp rise immediately prior to perias- 



3D hydrodynamical modelling of r] Carinae 



17 



3.5e-10 



3e-10 



2.5e-10 



1.5e-10 



5e-11 



^ 5*- • i .... / 

' A.*; *'"- '-^^ . ■ i- . J V " 

M ^ ^ ! ^ ^ 



Orbit-IA ' 
Orbit-RD 
RXTEdata 



hi 

I 



1 



1.2 



1.4 



1.6 



1.8 2 

Orbital Phase 



2.2 



2.4 



2.6 



2.8 



3.2 



3e-10 
2.5e-10 

2e-10 
1.5e-10 

1e-10 

5e-11 






Orbit-IA 




Orbit-RD 




; RXTE cycle 1 















3e-10 
2.5e-10 

2e-10 
1.5e-10 







Orbit-IA 






Orbit-RD 






/ 1 RXTE cycle 2 






/ ^ 1 


















! :', 












V 1.: 









3e-10 
2.5e-10 

2e-10 
1.5e-10 



Orbit-IA ■■ 
Orbit-RD 
RXTE cycle 3 ■■ 



0.9 



1.05 1.1 



Fig. 16. — 2-10 keV X-ray lightcurves for models Orbit-IA and Orbit-RD plotted against the RXTE lightcurve over three cycles (top 
row) and against the individual minima (bottom row - left to right: 1998, 2003.5, and 200 9). The estimated cosmic background has been 
removed from the RXTE data (^Corcoran et al...200ll : ICorcoran|[2005l : ICorcoran et al.llMIOh . 



2-5 keV ■ 7-10 keV ■ 2-10 keV 



Fig. 17. — Ray-traced images showing X-ray emission (erg s ^ cm ^ keV ^ ster) from model Orbit-RD at ^ = 1.0 in the 2-5 (left), 7-10 
(middle), and 2-10 keV (right) energy bands. The centre of the image is aligned with the centre of the simulation box and the viewing 
angles are i = 42° and 6 = 20° (§ 15. 2|) . The plots show a region of ±4 x 10-*^^ cm - large axis tick marks correspond to a distance of 
1 X 10^^^ cm. 



tron which is followed by a dip, and then there is a sec- 
ond smaller peak. As shown bv lParkin et alJ ()2QQ9f ). the 
first peak {(j) = 2.00) corresponds to a rapid increase 
in absorption as the emission region is obscured by the 
dense unshocked primary's wind, and the second peak 
{(j) = 2.02) corresponds to sightlines to regions of high 
intrinsic luminosity becoming closely tangential with the 
leading arm of the WCR (see the (j) = 1.02 snapshots in 
Figs. [T2l and [T3|) . The differences between the time at 



which the EWC in the models and the data begin to in- 
crease as periastron is approached, and also the time of 
the peak in column, suggest that our adopted azimuthal 
angle ^ = 20° is too small - if the value were increased 
slightly then the peak in column density would occur at 
a later orbital phase. 

We note that in the model of Oka zaki et al.l (|2008f ), 
which assumed the emission to originate from a point 
source situated close to the apex of the WCR, it was 
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Fig. 18. — Snapshot of the orbital (x — y) plane from model 
Orbit-RD at </> = 1.0 showing the 2-5 keV X-ray emissivity 
(erg s"-*^ cm~^ keV"-*^ ster). The region of postshock gas in the 
trailing arm of the WCR responsible for the excess of 2-5 keV X-ray 
emission is labelled "A" . The labels "B" and "C" refer to suggested 
sites for S and Si X-ray line emission. The position of the primary 
and companion star are indicated by the "PS" and "CS" labels, 
respectively. The arrow in the lower left of the image indicates the 
direction vector (in the orbital plane) for our adopted line of sight. 
The plots show a region of zbl x 10-*^^ cm - large axis tick marks 
correspond to a distance of 1 x 10^"^ cm. For the corresponding 
density, temperature, and speed images see Fig. 1131 

this secondary peak in column density at ^ = 2.02 which 
allowed the extended X-ray minimum to be reproducecfl. 
This is evident if one compares the position of the apex 
of the WCR relative to the leading arm of the WCR in 
the (j) = 1.02 and 1.04 snapshots in Figs. [T2l and [T3l 

5.2.4. Spectra 

X-ray spectra of r] Car are particularly useful as they 
put tight constraints on the energy dependence of the 
X-ray emission from the postshock gas, which constrains 
the postshock gas temperature and the preshock velocity. 
Therefore, the spectra provide valuable clues as to the 
variation of the preshock companion's wind speed as a 
function of orbital phase. 

In Fig. [20] we compare the X-ray spectra from mod- 
els Orbit-IA and Orbit-RD to XMM-Newton s pectra at- 
tained around periastron passage by Hama guchi et al] 
(|2007[ ). To facilitate our comparison the non- variable 
emission components from the outer-ejecta, the X-ray 
Homunculus nebula, an d the central cons tant emission 
component identified by iHamaguchi et all have been re- 
moved, so that the remaining emission is almost en- 
tirely contributed by the WCR. At ^ = 1.924 the 
match between models Orbit-IA and Orbit-RD and the 
XMM-Newton spectrum appears quite good with only 
a slight overestimation (underestimation) of the flux at 
E < 2 keV {E 3 — 7 keV). The relative normalization 
of the models and the XMM-Newton spectra suggests 
that our adopted companion wind mass-loss rate is rea- 

^ There was an additional source of absorption from a pile-up of 
unshocked primary wind gas in the SPH model of Okazaki et ak] 
(|20081 ') . In our grid-based hydrodynamical models we do not see 
such a build-up as this gas is processed through the shocks and 
thereafter flows downstream within the WCR. 



sonable, whereas the slope of the 7-10 keV flux indicates 
that a slightly lower preshock companion wind speed is 
required at this phase. 

The XMM-Newton spectra at = 1.988 and 1.990 lie 
at the bottom and top of a flare-like feature, respectively. 
Therefore, the fact that both model Orbit-IA and Orbit- 
RD overestimate the flux level at (/) = 1.988 is unsurpris- 
ing. Ideally, we would then like the models to slightly 
underpredict the observed flux level at (/) = 1.990 so that 
the average flux is approximately correct (i.e. neglecting 
the mechanism that produces the flare-like variability in 
the lightcurve). However, as the models are unable to re- 
produce the fiare-like features in the observed lightcurve 
(Fig. [T6|) . we take the level of agreement as reasonably 
good within the scope of the current investigation, and 
defer a detailed study of the flare-like features to future 
work. 

Comparing the goodness of fit between the XMM- 
Newton spectra an d the models presen ted in this work 
with the results of iParkin et al.l ()2009l ) we see a num- 
ber of interesting details. Firstly, at orbital phases prior 
to periastron in cycle 2 (0 = 1.924,1.988, and 1.990) 
IParkin et al.l showed that a reduced preshock compan- 
ion star wind speed improved the agreement between 
the observed spectra and their model. However, we now 
find that even without any inhibition of the preshock 
companion's wind (model Orbit-IA) the fits to the spec- 
tra are considerably better. This implies that a more 
accurate description of the hydrodynamic structure of 
the WCR around periastron (i.e. the highly asymmet- 
ric temperature distribution in the leading and trailing 
arms - see Figs. [12] and [T3|) helps to account for the be- 
haviour of the observed X-ray emission. In addition, our 
new model with modest radiative inhibition matches the 
spectra about as well as model Orbit-IA. This is because 
prior to (j) = 0.990 radiative inhibition of the preshock 
companion's wind is actually relatively minor. 

The spectra taken at orbital phases during the ex- 
tended X-ray minimum of cycle 2 (^ = 2.009, 2.015, 
2.018, and 2.023) are not well matched by the models, 
which overpredict the observed 7-10 keV flux by roughly 
an order of magnitude. Also apparent is the fact that 
the spectra do not drop off sharply as one tends to- 
wards lower energies. Examining the broadband images 
at these orbital phases shows that a significant amount 
of soft {E < 3 keV) X-ray emission originates from the 
tail of the WCR (see e.g. Fig. [TTj). In essence, there are 
two spatially discrete emission regions: the postshock gas 
near the apex of the WCR from which we only observe 
a significant flux of hard {E > 5 keV) X-rays due to 
the obscuration by the primary's dense wind, and the 
downstream gas in the tail of the WCR which lies in the 
wake behind the companion star. The additional softer 
X-rays in model Orbit-IA at = 2.009, 2.015, 2.018, and 
2.023 compared to model Orbit-RD is therefore due to 
the higher level of instabilities in the tail of the WCR 
in the case of the former. Such instabilities perturb the 
dense layer of postshock primary wind into the path of 
the postshock companion's wind and in so doing they 
help to reheat the gas through a series of shocks - this is 
well illustrated by a comparison of the = 1.020 snap- 
s hots in Figs. [121 and IT31 

IParkin et al.l ()2009[ ) found that an improved agree- 
ment between the model and data could be attained 



3D hydrodynamical modelling of rj Carinae 



19 



E 

° 1e+23 



Orbit- 1 A ■ 
Orbit-RD ■ 



H07-2-10keV [^— . i\ 



□ 



□ 



a 



1.4 1.6 
Orbital Phase 



Fig. 19. — Emission weighted column density (EWC) for models Orbit- IA and Orbit-RD for or bital phase range (j) = 0.9 — 2.1 (left panel) 
and 1.9-2.1 (right panel). For comparison, the column densities derived bv Ham aguchi et al.l (|2007l ) for the 2-10 keV X-ray emission are shown 
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at ^ = 2.009 by either strongly inhibiting the preshock 
companion's wind to a velocity of 1600 km s~^ (not 
tested in this work), or via a disruption/collapse of the 
WCR. The net effect of these alterations was to reduce 
the 7-10 keV flux by a factor of roughly 10. Interest- 
ingly, a collapse/disruption of the shock apex may also 
help to resolve the problem with the excess emission at 
E < 3 keV by later modifying the conditions in the WCR 
further downstream. We note again that although a dis- 
ruption/collapse of the WCR does not occur in model 
Orbit-RD, our adopted formalism for the coupling of the 
radiation fields to the opposing wind (Eq. [T2]) provides 
a lower limit to the degree of radiative inhibition of the 
companion's wind by the primary's radiation field. For 
instance, if we were instead to adopt k2 and a2 for the 
coupling of the companion's wind and the primary's ra- 
diation field the coupling would be stronger, and likewise 
the decelerative force would be greater. Therefore, the 
preshock companion's wi nd would be more s trongly in- 
hibited (see figure 21 of iParkin et al.l l2009l ) and post- 
shock gas would cool more rapidly making conditions 
more favourable for instability growth, and thus a dis- 
r uption/collapse of the WCR. 

iCorcoran et al.l (|2010f ) compared Chandra X-ray spec- 
tra taken at roughly the same point during the 2003.5 
and 2009 periastron passages, i.e. (j) = 2.03 and 3.03. 
Interestingly, the (j) = 3.03 spectrum showed an approx- 
imately identical shape to that at (f) = 2.03, albeit with 
a factor of 12 higher flux level. Encouragingly, the 
models agree much better with the (j) = 3.03 spectra 
and support our conjecture that the lack of an extended 
WCR disruption/collapse is consistent with the recent 
periastron passage. 

6. DISCUSSION 

The simulations presented in § H] of static stars at 
a separation corresponding to periastron (model Peri- 
RD) show that if the companion wind is suppressed suf- 
ficiently then the WCR displays catastrophic instabili- 
ties where dense clumps may hit the companion star. 
However, when orbital motion is included (§ [5l model 
Orbit-RD) the reduction in the preshock companion 
wind speed caused by radiative inhibition is smaller and 



higher preshock velocities stabilize the WCR. If one com- 
pares the results of simulations Peri-RD and Orbit-RD 
there appears to be a smoking gun for explaining the 
observed extended X-ray minima. The key lies in the 
radiative cooling of the postshock companion wind. We 
have argued that this can be brought about by the sup- 
pression of the preshock companion wind by radiative 
inhibition which, given our adopted relatively weak cou- 
pling between the primary's radiation field and the com- 
panion's wind, may not actually be such a hard thing 
to achieve. For example, the coupling between the radi- 
ation fields and the winds which we have adopted es- 
sentially provides a lower limit to the level of radia- 
tive inhibition which will be felt by the companion's 
wind (see § [3] and Fig. [3]). In reality, the couphng be- 
tween the radiation fields and the winds may lie some- 
where between the two extremes of each star, and in that 
case a greater inhibiting force would be imparted on the 
presh ock companion's w ind by the primary's radiation 
field (jParkin et al.l 120091 ). In addition, observations of 
massive star winds show them to be inhomo geneous (e.g. 
Fullerton et al. 1996|; iMarkova et all 120051 : iDavies et all 
20051 : iMoffat 2008 ). and the processing of clumps in the 
WCR could provide a sufficient density enhancement to 
trigger radiative cooling of the postshock companion's 
wind. Therefore, a useful direction for future investiga- 
tions would be to perform simulations focused on perias- 
tron passage which explore a broader range of radiation- 
wind couplings and also the influence of inhomogeneous 
winds on the postshock gas conditions. 

The factors which we mention in relation to trigger- 
ing a disruption/collapse of the WCR may also help to 
explain the observed fiare-like rapid va riations in X-ray 
brightness. iMoffat fc CorcoranI ()2009l ) presented sepa- 
rate models based on clumps from the primary's wind en- 
tering the WCR and on an unstable WCR and found that 
each model had its fair share of difficulties reproducing 
the observations. Based on the results of models Peri-IA 
and Peri-RD, i.e. oscillations in the WCR causing rapid 
variations in the X-ray flux (Fig. [6]), re-examining the 
combined role of instabilities and wind clumping using 
hydrodynamical models may be worthwhile. Further- 
more, although relatively small clumps can be efficiently 
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destroyed in largely adiabatic WCRs (|Pittardll2007[ ), this 
will not be the case for a radiative WCR. Large clumps 
traversing a highly unstable WCR consisting of radia- 
tive shocks may cause much greater variability, and could 
therefore account for the onset of large amplitude flare- 
like features as periastron is approached. 

In model Peri-RD the clumps which collide with the 
companion star are removed from the computational do- 
main when the wind is initiated after every time step. 
In reality, these clumps might disrupt the generation 



of the wind. While aspects of thi s scenario a re rem- 
iniscent of the model proposed by iSokerl (120051) . there 
are some crucial differences. In the ISokerl model the 
clumps originate in the preshock primary wind, whereas 
in model Peri-RD the clumps are formed by the frag- 
mentation of the dense layer of postshock primary wind 
by the NTSI. Furthermore, in our work clumps do not 
traverse the shocks due to gravitational attraction, but 
instead are catapulted by violent instabilities. Since 
in our simulations the preshock stellar winds are mod- 
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elled as homogeneous, w e cannot directly test "SokePs 
model (althou gh see Wald er fc Foli ni 2002; Pittard 2007| 
iParkin fc Pitt ard 2010, for a discussion of the ablation 
of clumps in postshock flow). 

Ideally, one would like to resolve the region around the 
stars such that the critical point radius of the wind accel- 
eration region is well sampled. It is, however, a compu- 
tationally demanding task to resolve the highly eccentric 
orbit of the stars and also the critical point radius (for 
the companion star this was ^ 1.03 i^*2)- It may there- 
fore be advantageous in future models to adopt a sub-grid 
treatment of wind accretion which modifies the driving of 
the winds. For instance, "switching off" the initiation of 
the wind over those parts of the stellar surface in which 
the ram pressure of incident clumps exceeds that of the 
wind. 

Photo-ablation by the stellar radiation fields clearly 
has important consequences for the fate of cool (~ 
10^ K), dense postshock primary wind. With the stars 
in T] Car being so immensely luminous, photo-ablation 
see ms ine vitable. However, in this work we adopted 
the ICAKI formalism for radiative driving, whereby the 
radiation- wind coupling is calibrated by the stars ability 
to drive its own wind. Therefore, the quantitative ac- 
curacy of future models would benefit from the use of a 
more suitable treatment of photo-ablation of circumstel- 
lar gas. 

7. CONCLUSIONS 

Three dimensional hydro dynamical simulations of the 
enigmatic super-massive binary star system 77 Car have 
been presented which include the radiative driving of the 
stellar winds, radiative cooling, gravity, and orbital mo- 
tion. A suite of simulations were then performed to ex- 
plore the role of wind acceleration, interacting radiation 
fields, and instabilities in the WCR on the gas dynam- 
ics and resulting X-ray emission. We summarize our key 
conclusions as follows: 

• When one separates the coupling between the ra- 
diation fields of the stars and the opposing wind 
the level of radiative inhibition of the compan- 
ion's wind is reduced. Consequently, in static- 
star calculations with a periastron separation the 
companion's wind attains a preshock velocity of 

2200 km s~^ - notably higher than the estimate 
of 1500 km s~^ when the radiation- win d cou- 
plings are not separated (jParkin et al.ll20Q9l ). 

• Despite a reduced level of radiative inhibition, 
when the wind acceleration regions are considered 
(in static-star calculations with a periastron sepa- 
ration) the WCR becomes massively disrupted by 
non- linear thin shell instabilities (NTSIs). This oc- 
curs because the lower preshock companion wind 
speed acquired due to wind acceleration and ra- 
diative inhibition leads to a lower postshock gas 
temperature. As such, radiative cooling becomes 
important in the postshock companion's wind and 
thermal pressure - which prevents the growth of 
the NTSI - is lost. The WCR is distorted to such 
an extent that dense fragments of cool postshock 
primary wind are repeatedly driven deep into the 
companion's wind acceleration region and in some 
instances collide against the star. 



• When orbital motion is included the catastrophic 
disruption of the WCR is not reproduced. The 
root cause of this difference lies in the rapid orbital 
motion of the stars around periastron which acts to 
increase the preshock velocity, and thus postshock 
pressure, of the companion's wind. In so doing the 
stability of the WCR against thin-shell oscillations 
is increased. The influence of orbital motion will 
depend on the adopted system parameters. 

• Large-scale, high resolution simulations show a 
number of interesting dynamical effects, including 
the influence of the stellar radiation fields on the 
growth of instabilities at the shocks, the photo- 
ablation of the cold dense layer of postshock pri- 
mary wind, and the highly asymmetric tempera- 
ture distribution in the arms of the WCR. 

• The models provide a reasonable match to the ma- 
jority of the RXTE lightcurve and, compared to 
previous models, the X-ray spectra agree much 
better with XMM-Newton observations obtained 
just prior to periastron. However, the extended 
X-ray minima are not reproduced by model Orbit- 
RD (which uses a relatively weak coupling between 
the primary's radiation field and the companion's 
wind) and the 7-10 keV X-ray emission is overes- 
timated by roughly an order of magnitude. Yet, 
when the shock is heavily disrupted (in the static- 
stars simulation), the 7-10 keV X-ray luminosity of 
the postshock gas is an order of magnitude lower 
than the undisrupted case. This shows that dy- 
namical instabilities in the WCR could be a key 
mechanism in explaining the X-ray observations. 

• From a comparison between the model and XMM- 
Newton column densities around the 2003.5 perias- 
tron passage, a reduction of the stellar wind mass- 
loss rates by a factor of ~ 2 is implied. The re- 
vised mass-loss rates are 2.4 x lO~^M0yr~^ 
and (7 — 12) x 10~^ yr~^ for the primary and 
companion wind, respectively. However, these es- 
timates do not account for the possible presence of 
highly ionized, and thus lower opacity, gas close to 
the stars, and seem unlikely given other evidence 
for m ass- loss rate s grea t er than those assumed 
here (iHillier et ap I2001L 120061 : Ivan Boekel et al.l 
120031 : iGroh et al.ll2010al ). Alternatively, it may be 
that the observationally derived column densities 
around periastron are subject to uncertainty due 
to fitting si mpler plasma models to complicated 
systems (see iPittard fc ParkinI I201Q[ ). and/or the 
process of disentangling different spectral compo- 
nents. 

• To reproduce the extended minima likely requires 
effective radiative cooling in the postshock com- 
panion's wind close to the apex of the WCR. This 
could be achieved by stronger radiat ive inhibition 
of th e preshock companion's wind (jParkin et al.l 
|2009^, albeit with a stronger radiation- wind cou- 
pling than explored in this work. If the compan- 
ion's wind is sufficiently suppressed, or oscillations 
in the WCR are especially vigorous, a catastrophic 
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disruption/collapse of the WCR onto the compan- 
ion star may occur. 



We close with a note that the differences between the 
recent X-ray minimum and those previously observed 
could be accounted for by the stochastic nature of the 
growth of non-linear instabilities in the WCR. For ex- 
ample, we can speculate that the growth of instabilities 
in the 1998 and 2003 minima sufficiently disrupted the 
WCR to result in an extended quenching of the 7-10 keV 
X-ray emission. In contrast, during the shorter 2009 
minimum the WCR was contorted by instabilities to a 
lesser degree, perhaps due to less significant radiative 
cooling of the postshock companion's wind, so that ul- 
timately an extended disruption and/or collapse of the 
WCR against the companion star did not occur. This 
hue of reasoning suggests that the differences between 
the most recent and the previous X-ray minima may be 
divided by a fine line related to the postshock gas con- 
ditions, and is an intriguing possibility given that mul- 



tiwavelength observations suggest recent changes in the 
wind of 77 Car (Martin et al. 2010: Mehner et al.ll2010bl : 
iCorcoran et al.lfeoiOl) . 
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APPENDIX 

A. INITIATING THE INSTANTANEOUSLY ACCELERATED WINDS 

In simulations where the winds are instaneously accelerated we do not need to resolve the wind acceleration region. 
Therefore, the only resolution requirement comes from having a sufficent number of cells between the stars to accurately 
model the WCR dynamics. With this in mind, computational resources can be saved by making the radial distance 
into which the winds are initiated i^map a function of the separation of the stars, 

i^map = i^peri X 2^-^^"^^ (Al) 

where i^peri is a radial distance corresponding to approximately 10% of the distance between the stars at periastron 
(> 6 cells), 1(f) is the number of nested levels of refinement required at a given orbital phase to accurately model the 
fiow dynamics between the stars, and /pen = /((/> = 0.). 

B. REFINEMENT CONDITION FOR POSTSHOCK GAS 

To ensure that postshock gas is sufficiently well resolved to accurately describe the WCR dynamics we have im- 
plemented an additional criterion for grid refinement. Given the desired number of cells between the stars, ngep, the 
number of nested levels of refinement required at a given orbital phase 

h = log2 (^) , (Bl) 

V ^base / 

where Creq = (isep(^)/^sep IS the resolution required, (isep(</^) is the separation of the stars, and Cbase is the coarse grid 
resolution. For models Orbit-RD and Orbit-IA, ngep — 200 cells. 

C. ESTIMATING THE INFLUENCE OF ORBITAL ACCELERATION ON THE WIND 

The motion of a star in an eccentric orbit introduces centripetal, ^cen, and line-of-centres, ^loc, accelerations. Con- 
sidering the fiow along the line of centres, 

and, 

^orb 27raecosa; 

where Vorh is the orbital velocity, R is the radial distance to the system centre of mass (= a(l — e^)/(l + ecoso;)), 
uj is the true anomaly (= 27r(/)), Morb is the mass term relevant to the particular orbit under consideration (e.g. for 
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Fig. 21. — The orbital acceleration of the companion star around periastron passage (left panel) and the radiative line force calculated 
for the companion's single star wind profile (right panel). The components for the centripetal and line-of-centres, as well as the total are 
shown in the left panel. Note that these values are calculated in the frame of reference of the companion star (c.f. Fig |14|) . 
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Fig. 22. — Companion star wind velocity along the line of centres between the stars. The primary star is situated at a distance of 
359 R0 (corresponding to periastron separation when e = 0.9). The following solutions are shown: without radiative inhibition and orbital 
motion (No RI)(i.e. single star), without Rl and with orbital motion (Orbit + No Rl), with Rl and without orbital motion (Rl) , an d with 
Rl and orbital motion (Orbit + Rl). Dynamically consistent FDCFs have been used in these calculations (c.f. Eos. ICSl and [C9| . The 
orbital parameters are noted in Tabled and stellar and line driving parameters are noted in Tabled 



the barycentric orbit of the companion star Morb = GMf / [Mi + M2)^), and P is the orbital period. The left panel 
of Fig. [2T] shows ^cen and ^loc evaluated for the companion star of rj Car. Comparing the acceleration due to orbital 
motion against the radiative line force calculated for the companion star one sees that the latter has a higher magnitude 
close to the star where the mass- loss rate is set (right panel of Fig. [2T]) . 

The influence of the orbital acceleration on the mass-loss rate and velocity profile of the wind can be estimated by 
incorporating ^cen and ^loc into the equation of motion 



F(r, dv/dz) = 



dv 
dr 



dAlj_ 

dr 



2s^ 

r 



9ra.d 9orh') 



(03) 



where r is the distance from the centre of the star along the line-of-centres, s is the speed of sound, tjj is the gravitational 
potential, and ^orb = ^cen + 5'ioc is the orbital acceleration. The radiative line force, ^rad takes the form of Eq. [TTl 
The value of Vorh at a given radius will depend on the time at which that parcel of gas left the star. However, as 
we are mainly interested in the affect of orbital motion on the base of the wind (i.e. where we initiate the winds 
in the hydrodynamic simulations), we neglect the radial dependence of '^orb, and take '^orb( ^) ^ ^'or b • A s the orbital 
acceleration is applied to the base of the wind we replace the R in the denominator of Eqs. ICll and [C2l with R — R^. 
Note that we do not use R — r in the denominator as this would introduce an unrealistic singularity at the c entre of 
mass of the binary system. To proceed, we make the coordinate transform using the substitution of variables (jAbbottI 

[Hsoi), 



-2GMi(l-ri) 



dr 



(04) 
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leading to 



where 



F{u^ w') 



- U;' + h(u) - CB(u, w, w')w''', 
w ' 



K2{u,w,w')A{u), 



and 



A{u) = 

with C given by equation (13) of [Stevens Pollock! ([T 



Ki{u^ w') 



U — Ud^ 

. The FDCFs 

-j^ f _u_\ 2w f _u_ 

J w'u \ 



(C5) 



(C6) 



(C7) 



1+a 



and 



K2{u,w,w') 



(C8) 



1 + Q! 



(C9) 



where Ud 
function, 



with 



'^(<^sep), = u{R^). The orbital motion terms, which only depend on r, are contained in the potential 



h{u) = 1 + 



M2(i-r2) 
Mi(i-ri) 



A{u) 



1 ^orb^* 



^orb 



■C, 



2 2 / 27rae'yorb cos uj o 



(CIO) 



(Cll) 



where i^orb = u{R). The critical point of the wind profile is solved for using: i) the equation of motion F{u^ w') = 0, 
ii) the singularity condition dF/dw' = 0, and iii) the regularity condition dF/du + w'dF/dw = 0. The velocity 
gradient terms present in the FDCFs (Eqs. [C8l and IC9[) require the solution for the wind profile to be attained using 
an iterative process, whereby mon otonic FDCFs are used to attain an initial solution, then velocity gradient terms are 
included in subsequent iterations (jPauldrach et al J 119861 ). 

Fig. [22] shows wind solutions for 77 Car with/without radiative inhibition and orbital motion, and calculated with a 
stellar separation corresponding to periastron. In these calculations ^rad has the form of Eq. [TTl whereas in the hydro- 
dynamic simulations ^rad has the form of Eq. [121 Hence, these calculations are intended to illustrate the qualitative 
effect of orbital acceleration on the wind. The inclusion of orbital motion slightly steepens the velocity profile close 
to the star. This is the opposite effect to stellar rotation, which causes a shallower velocity profile (Abbott 19^. 
Interestingly, neglecting radiative inhibition, the terminal wind velocity is lower when orbital motion is included. How- 
ever, when radiative inhibition is included, and recalling that the WCR occurs at 100 R© from the companion star 
(at periastron and along the line of centres) in model Orbit-RD, there is a minor increase in the preshock velocity. 
Examining the critical point radius and mass- loss rate we find reduction factors of 0.12 % and 2.3 %, respectively. 
Therefore, orbital acceleration contributes a negligible difference in comparison to the interplay between the radiation 
fields and the gravity of the stars. 



